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Abstract 

We compute the renormalization mismatch displayed in 1-loop approximation 
by classically equivalent 4-quark operators and coming from different possible 
definitions of the 75 matrix in dimensional regularization. The result is then 
employed to study the effect of the various treatments of 75 upon the size 
of radiative corrections to 4-quark condensates in the QCD sum rules for p 
and A\ mesons. We find that a fully anticommuting 75 which automatically 
respects non-anomalous chiral Ward-Slavnov identities leads to considerably 
smaller corrections and reduces theoretical uncertainty in the QCD prediction 
for the r hadronic decay rate. 



1 Introduction 



Dimensional regularization (DR) |TJ coupled with the minimal subtraction scheme (MS) 
is widely acknowledged as theoretically sound and extremely valuable in practice calcu- 
lational framework. In a theory like QCD with parity conserving amplitudes it respects all 
the symmetries of the theory such as gauge invariance, Bose symmetry and Ward-Slavnov 
identities. In addition, dimensionaly regularized divergent Feynman integrals are in many 
respects very similar to the convergent ones, which allowed to develop quite powerful 
methods of their analytical calculation. 

A notorious shortcoming of DR is the absence of a natural generalization of the Dirac 
matrix 75 to noninteger values D 7^ 4 of the dimension of space-time. Indeed, the well- 
known relation 

t r (75^ '[111 'n2^l '^/) Th) — _ 4^ e /il/J2M3/i4 (1) 

which in fact fixes 75 in 4 dimensions, is not compatible with another property of 75 
valid at D = 4 

{75,7m) = 0, (2) 

The problem is that the totally antisymmetric e-tensor is a purely 4-dimensional object 
and thus can not be self-consistently continued to D dimensions. It may be shown that if 
(0) is to be held literally for all values of fi = 1 . . . D (and the standard properties of the 
trace operation like cyclycity is respected) then the rhs of ([]]) should vanish identically 
||. Such a defined 75 does not, hence, go smoothly at D — > 4 to the 4-dimensional 

# = -27V7V (3) 

in contradiction with general requirements to which every self-consistent regularization 
must meet [|J]. From the same point of view the relation (|2]) may be freely violated for 
the values of index \i lying outside of the physical 4-dimensional subspace. 

The existing solutions of the problem may be divided into two classes. The first one 
favours a nonvanishing anticommutator {75, 7^} at fi ^ 0, 1, 2, 3 and is discussed in general 
in H and || . An explicit version of such a scheme was proposed by 't Hooft and Veltman 
and elaborated by Breitenlohner and Maison |], |3|] . Another class keeps (§) intact at least 
for amplitudes without traces with odd numbers of 75 matrices while trying to reproduce 
the identity ([3]) either by hand J?| or, in a bit more sophisticated manner, by redefining 
the trace operation for the traces with odd number of 75 || [|. 

Leaving aside a (rather academic) task of finding an explicit formal proof of self- 
consistency of these solutions, we will study in this work another aspect of the 75 prob- 
lem — different possibilities of choosing between classically equivalent 4-quark operators 
which begin to differ in higher orders due to different treatment of the 75 matrix. As we 
shall see below, different choices lead to different theoretical predictions of the QCD sum 



rule method |T(J because of a commonly accepted approximate procedure of the vacuum 
saturation in estimating the vacuum expectation values of 4-quark operators. On consid- 
ering a number of particular physical problems, it will be demonstrated that the choice of 
fully anticommuting 75 considerably reduces theoretical uncertainty coming from higher 
order corrections. 
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2 The 75 ambiguity in 4-quark operators 

The problem under discussion appears even in cases when initial amplitudes comprise no 
75 matrix at all. Indeed, let us consider the standard derivation of the QCD sum rules 
for the p and A\ mesons |10|. It starts from constructing the Wilson expansion [11| for 
the correlator 

IW<z) = i J dxe^ < T 3ll (x)jt(0) > = {q,q v - g^q 2 )Tl(Q 2 ) (4) 

n(Q 2 ) = c + Ci >0 + cf ' A ^ ° l >0 + ■ . . (5) 

where = u 7^(75)0? is the current for the p(Ai) meson. 

The power corrections of order 0(1/Q 6 ) to the rhs of the expansion ([5]) come from 
4-quark condensates (that is the vacuum expectation values (VEV) of 4-quark operators) 
and lead to most important numerically non-perturbative corrections to the sum rule for 
the p meson. In the lowest order of perturbation theory the contributions of the 4-quark 
operators are well-known and read [TIJ 

= -\{u lfl t a u + d^t a d) ■ (* 7 W) - 2u m^dd^^u (6) 
9 " 

for the vector correlator and 

C l 6 ' A Ql 2 



9 2 9 



u ^t a u + d 7 M t a d) • Yt a, &) - 2u^t a d d^t a u (7) 



for the axial vector one. Here t a , a = 1, ... 8 are the colour Gell-Mann matrixes normalized 
as Tr(t a t b ) = ^5 ab ; (^fT^f) is a shorthand for Y^f= u ds(^f^^f)] 9 * s the quark gluon 
coupling constant such that a s = g 2 /An. 

Let us consider in detail the structure of the second term in (||) proportional to the 
4-quark operator 

0[=u l5llM t a dd l5llM t a u. (8) 
In fact, the very calculation produces a bit different operator 

1 = uT^t a ddT^t a u, (7^7,-7,7^) (9) 

while the result @ is obtained after a use of the identities 

r [3] ( g )r [3] = 6 75 7m(S)757m- (10) 
Similarly, the initial form of the second term in (|7|) 

0' 2 = u-i^t a dd-i^t a u. (11) 
is (with an appropriately changed coefficient function) 

Oi = ur® l5 t a ddrW l5 t a u, (12) 
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The relation (|T0| ) is valid in 4- dimensional space due to the very definition of the 
matrix (||). However, it loses any operational sence at generic within the framework 

of dimensional regularization. This implies that if both operators 0\ and 0[ are to be 
minimally renormalized, then the difference 

5Q 1 = O x - QO[ = 0{a s ) (13) 

is an evanescent operator [tH] which vanishes in the tree approximation but typically 



receives nonzero contributions in higher orders. This also means that higher order correc- 
tions to the coefficient functions (CF) Cg also depend on specifying which operators are 
chosen as basic ones. It should be stressed that this 75 ambiguity comes from the absence 
of a unique canonical representation of 4-quark operators in dimensional regularization, 
which, in turn, does not allow to define the idea of minimal subtraction of UV poles in 
an unambiguous way. 

The above discussion does not imply in any means that the theoretical predictions for 
the contributions to the rhs of ([|) from 4-quark operators are ill-defined in principle: the 
VEV's of the latter should, of course, also depend on the choice of the renormalization 
prescription in such a way that the whole combination J2i C 1 qO\ is invariant as it must be. 

The real problem lies in our present inability to find the VEV's of 4-quark operators in 
a sufficiently accurate way. Indeed, the common practice is to use the vacuum saturation 



(VS) procedure fLOf. It states that the VEV of a 4-quark operator O = (^Ti^)(^T 2 ^] 



may be estimated according to the following formula 

< ° >0 = lii t Tr ( r i) Tr ( r 2) - Trpfo)] -<qq> 2 (14) 

where qq = uu, dd or ss. 

It is clear that this procedure is too rough and does not feel in any way such fine details 
as the concrete mode of the renormalization of a 4-quark operator. Hence, if one would 
like to estimate the effect of higher order corrections to the CF's C l & the numerical results 
will do depend on specifying 4-quark operators in D dimensions. Unfortunately, these 
corrections prove to be of considerable size (see \Y1\ and below) at least in the operator 



basis featuring objects like r' 3 ' (^r' 3 l Let us try to search for a distinguished operator 
basis for which the VS would presumably work better. 

Of course, from a general point of view there is an infinitely large variety of possible 
choices of a (different at D 7^ 4 and equivalent at D = 4) 4-quark operator. For instance, 
one could mix operators 0\ and 0[ as follows 

0" x = O x ■ sin 2 a + 60[- cos 2 a (15) 

with an arbitrary angle a G {— 7r/2, 7r/2}. 

On the other hand a lot of calculational experience invariably demonstrates that if 
amplitudes under investigation do not comprise traces with an odd number of 75 at 
all then the use of a fully anticommuting 75 is quite self-consistent and leads to great 
simplifications in doing practical calculations. Much more important, this choice is in 
a sense unique for the amplitudes from above described class (these will be referred to 
non-singlet ones in the following) as it respects the chiral SUjj(3) ® SUr(3) symmetry of 
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the QCD lagrangian. This is the case since on the purely diagrammatic level the chiral 
symmetry relies on the possibility of freely anticommutating of 75 with 7^ at generic 
value of the index \i. It seems very reasonable to employ identities like (|IUD in order to 
maximally simplify kinematical structures appearing. Thus, we will treat the operators 
0[ and 0' 2 with completely anticommutating 75 as "natural" ones and reexpress through 
them the operators^ 0\ and O2. 
A simple calculation gives: 

AO a 471 3 2 a 471 ( 16 ) 
+- -r^u la dd la u + - -±{u iXu + d 7 a t a rf) (* T W), 
3 Att 9 An 

2 = (l-^)-60' 2 + 2 ^u lal5 t a d d 7 Q 75^ 
40 a 4?r 3 2 4 a ( 17 ) 

+- -r^u w 5 dd lal5 u + - -a-iu i a t a u + d rt a d) (^yv*). 

3 Air 9 Air 

To finish this section, a few technical words about the calculation. It is clear that 
the very difference Oi — 0\ vanishes in the classical limit and consequently may come 
only from UV divergences which manifest themselves as poles in e in DR |TJ|]. Thus, 
the calculation is conveniently performed hand in hand with the standard procedure of 



evaluating one- loop anomalous dimensions of 4-quark operators [TO (see also where 
similar calculations were performed in the context of evaluating QCD corrections to the 
effective weak Hamiltonian. 



3 QCD sum rules for the p and A\ mesons 

An explicit calculation of one-loop corrections to 4-quark condensates in @ gives |E2| (in 
the commonly accepted MS scheme L — In 4j ) 

1071 \ 

(u lfl t a u + d lfl t a d) ■ (* 7 W) 

u T [3 k a ddT [3] t a u 



Cj v Qi 
9 2 



1 

3 



1 + 



71 



7T 



95 r 107 

— L H 

72 48 

9 r 431 

7^L + 



96 



O's 



(16L-12) uT^ddT^u 
+ (30L-f) uT^ddY^u 



24vr 



(18) 
(19) 



(u 7 m 7 5 m + J 7 M 7 5 d)(* 7^75*) 

(u i^t a u + J 7 M 7 5 rd)(^ 7^75^*) 



• (20) 



■""Note that we are concerned only with 4-quark operators appearing in the tree approximation; for 
operators which show up hrst in 1-looop corrections the subtlety under discussion becomes relevant 
starting at two loops. 
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For the axial vector correlator the corrections look exactly the same, only the matrices 
r' 3 J and T M should be multiplied by 75 . After the use of the VS approximation we get 
for the p meson 

^iv 224vr / ,705 13 rS a s \ . 9 , 

Cf < Ol >„ = -_(! + (_+ -L)^) (a, < ft >>„) (21) 

and for the Ai meson 

a 3527T / .777 149 rX a»\ , 9 . 

At last, a straightforward application of (0)>(IIZD l ea ds us to a new result for flT8 - 20|) and 
consequently for that of (|2"T|): 

, iy 224vr / .685 13 _. a s 



Cr <Ol>, = — {! + <- + -L^ ) (a, < ft >^ (23) 



and for the Ai meson 



^f> a ,™ 3527T /_ .917 149 r ,a s \ 

Cl' A <0*> n = 1 + + L)— (a s < qq > 2 n ) 24 

e 60 81 V v 528 396 ; 7T ; v HH oJ 



Thus, we observe that uncomfortably large radiative corrections in (pi|)and ( P4 ) are trans- 
formed into quite moderate ones after the above described transition to the new basis of 
4-quark operators. Moreover, now the a s corrections get of approximately the same (rel- 
ative) size for the vector and axial correlators — a fact that will be of special importance 
in the next section. 

Note also that the terms proportional to L in ( fffi ) and (24) are quite small. To our 
opinion, this fact demonstrates that the effects due to two-loop anomalous dimensions of 
operators in ( p0|) are presumably small and may be neglected. On the other hand, the 
combination of operators a s {p) < qq > (fi) < qq > (//) has only very weak dependence 
on fx as its anomalous dimensions is fortunately rather small. This means that if the VS 
approximation is a good one at some scale fi, then it should be considered as a reasonable 
approximation at a diffrent scale fi'. 



4 The tau hadronic width 

As was pointed out on J16], |17], [Tj| some time ago, the methods of perturbative QCD can 
be applied to estimate the decay rate ratio 

R T = r(T r^ hadrons) . (25) 
r(r -> v T ev e ) 

An updated theoretical discussion of the ratio (plf ) within the QCD frameworks was 
recently given in [p0 |. 

Schematically, R T may represented in the following form (in the chiral limit of the 
massless u, d and s quarks and neglecting the logarithmic dependence of the Wilson 
coefficients) 

R T = R° {l + ^ + 5 .202(^) 2 + 26.37(^) 3 + S D = 6 + ...}. (26) 
[nn n J 
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Here the perturbative coefficients are appropriate for the number of flavours / = 3 and 
a s = a s (M T ) and the power correction of order 1/M® is directly related to the contribution 
of 4-quark operators to (|5[), viz. 

-£(<# VA1 )<0>i) 

6 D =* = {6p* + 5?*)/2, ^ = 24^—* w (27) 



As was argued in p0[ the only significant sources of uncertainty in the QCD prediction 
for R T are due to the uncalculated perturbative QCD correction of order aj(M T ) and due 
to the power correction of order 1/M®. In this section we will concentrate on the latter. 
The estimation of the power correction made in |20 amounts essentially to keeping 



the lowest order contribution to the CF 5 D 6 and using the VS approximation. The result 
is 

D=6 ( 7 \ 256tt 3 pa s < jj> >-' 
V ={-ll) "27 M| • (28) 

Here pa s < ipip > 2 m (3.8 ±2.0) x 10~ 4 GeV 6 is an effective scale invariant matrix element 



which is determined phenomenologically [21, 22, 23]. It is introduced in order to partially 



take into account deviations from the VS approximation. An important observation is that 
when the vector and axial vector terms are averaged to give S D=6 , there is a large mutual 
cancellation between these contributions. It is the assumption — that this cancellation 
also holds for the uncertainties of separate contributions^ — , which is crucial in producing 
a rather accurate prediction for a s given in f20|l, viz. 



a s (M T ) = 0.34 ±0.004. (29) 



This assumption is certainly quite optimistic and has already been criticized in p4| , [2J 
From formal point of view it is self-consistent provided one may neglect higher order 
corrections to the coefficient functions involved. However, the above discussion of the 
leading 1-loop corrections to the CF Cg V ^ A clearly demonstrates that this neglection is 
very questionable unless one deals with the primed basis of 4-quark operators and avoids 
objects like T' 3 ' ® r' 3 l in favour of those featuring the completely anticommuting 75. 



5 Conclusions 

In this work we have studied the problem of the dependence of radiative corrections to 
the coefficient functions of 4-quark condensates on the exact way of treating the 75 matrix 
and related objects specifying a 4-quark operator. It was argued that a choice of a fully 
anticommutating 75 matrix which respects the chiral symmetry of the QCD lagrangian 
is theoretically preferred and leads to better convergent perturbation series. It also gives 
an extra argument in favour of the stability of (]28f) and (p9| ) with respect to radiative 
corrections. 

We would like to thank J. H. Kiihn and A. Kwiatkowski for their support and helpful 
discussions. 

2 This would become true if the prescription ( ff^ ) could be made exact by multiplying its rhs with an 
(operator independent!) correction factor p. 
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